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We investigate the thermodynamic efficiency of sub-micro-scale heat engines operating under the 
conditions described by over-damped stochastic thermodynamics. We prove that at maximum power 
the efficiency obeys for constant isotropic mobility the universal law rj = 2 rjc/{4 — rjc) where rjc 
is the efficiency of an ideal Carnot cycle. The corresponding power optimizing protocol is specihed 
by the solution of an optimal mass transport problem. Such solution can be determined explicitly 
using well known Monge-Ampere-Kantorovich reconstruction algorithms. Furthermore, we show 
that the same law describes the efficiency of heat engines operating at maximum work over short 
time periods. Finally, we illustrate the straightforward extension of these results to cases when the 
mobility is anisotropic and temperature dependent. 
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Present-day technology makes possible high accuracy 
handling of colloidal micro and sub-micro particles to 
study experimentally thermodynamic processes out of 
equilibrium Q. A recent experiment, Q used a highly 
focused infra-red laser beam to confine a S/irn diame¬ 
ter colloidal particle into a parabolic potential of tunable 
stiffness. By alternately varying the stiffness of the po¬ 
tential and, on much faster time scales, the temperature 
of the solvent where the particle was suspended pro¬ 
vides a clear evidence of the possibility to construct a 
Stirling heat engine at the micro-scale. This experiment 
is ground-breaking because it opens the way of exploring 
the limits to which it is possible to scale down heat en¬ 
gines while retaining the same working principles of their 
macroscopic counter-parts. 

In order to set the scene for the present discussion it 
is expedient to start by recalling the protocol governing 
an ideal Stirling cycle. During such a cycle, a working 
system undergoes a sequence of four thermodynamic pro¬ 
cesses while being alternately in contact with two equi¬ 
librium reservoirs at different temperatures. First, an 
isothermal heat uptake Qm from the environment at 
temperature (measured in energy units) ■ During 
this process the system expands, i.e., it converts heat 
into work done on the environment. Second, an iso- 
choric (constant-volume) cooling of the environment to 
a temperature ■ Third, an isothermal com¬ 

pression during which the working system transfers an 
amount Qout of heat to the cold reservoir. Finally, an 
isochoric heating of the environment back to ■ The 
neat effect of the cycle is the extraction of an amount of 
work Wout from the system with efficiency measured by 
the ratio 77 = WoutlQin- General thermodynamic rea¬ 
soning 12| proves that the efficiency of any heat engine 
cyclically working between two classical equilibrium heat 
baths cannot exceed Sadi Carnot’s bound 


, Ph 
VC = k- ^ 


( 1 ) 


The bound is attained by performing a reversible cycle. 


Thus, Carnot’s bound can be classically attained in the 
same adiabatic limit when the power output tends to 
zero. 

Coming back to the experimental conditions imple¬ 
mented in [^, a distinctive property of micro and sub¬ 
micro scale heat engines is that they operate with char¬ 
acteristic energies of the order of thermal fluctuations 
ilH- The state Ct of working system at time t 
must be thus characterized by a probability distribution 
m(z, t)dz = P{z < Ct ^ ^ d- d^:) over the space D of ad¬ 
missible configurations. The expansion (compression) of 
the system corresponds then to the increase (decrease) 
of the fluctuations of the particle induced by varying the 
temperature of the solvent. Thermodynamic indicators 
such as heat and work also become fluctuating quanti¬ 
ties. The analysis of thermodynamic cycles needs to be 
rephrased in the language of stochastic thermodynamics 
(see e.g. [i^, the more recents [2^ and refs, therein). 

The micro-scale heat engine of relies on the heat 
cycle theoretical model 2^ to extract thermodynamic 
quantities from the measurements of the colloidal par¬ 
ticle position. In [2^ the dynamics of the micro-scale 
particle is modeled by Langevin-Smoluchowski (over¬ 
damped) equations. Quantitative analysis of the Stir¬ 
ling cycle is then accomplished in the special but impor¬ 
tant case of Gaussian fluctuations. The cycle starts at 
time when the system statistics is a normal distribu¬ 
tion centered at the origin and with assigned variance 
tj^. Isothermal expansion at temperature is enacted 
by requiring that at a further time t~ the system proba¬ 
bility be still normally distributed with zero average but 
with assigned variance cr^ > cr^. Then the probability 
remains continuous across the discontinuity of the tem¬ 
perature profile which at t = to drops to ■ Isothermal 
compression brings back the variance of the distribution 
to (Tq at t = tff. Finally, the cycle closes at t = t{ with 
an abrupt temperature transition from I3~^ to . In 
order to yield a finite power output the cycle operates in 
a finite period T = tf — it. Hence a natural definition of 




2 


optimal efficiency Q is that of the protocol maximizing 
the power output over the duration of the isothermal ex¬ 
pansions and the period T. The result for the Gaussian 
cycle is [ 2 ^ 


V* = 
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More qualitative considerations led [ 2 ^ to suggest for 
non-Gaussian statistics and anisotropic temperature- 
dependent mobility 


also [21| ) and subsequently applied to predict refinements 
to the second law of thermodynamics and to Landauer’s 
bound [ll which have found experimental confirmations 

3.g. 13|. 


see e.i 


Model. We suppose that the configuration space of 
working system is while its state evolves according 
to the Langevin-Smoluchowski dynamics 


dCt = dt + 



( 4 ) 


Vi. = 


VC 


2- avc 


( 3 ) 


with 0 < a < 1 admitting an explicit expression only 
in the case of isotropic mobility and recovering ([2]) for 
constant mobility. In order to compare these predictions 
with experiments, two questions were left open. First, to 
derive in the most general case the value of a in terms of 
the parameters of the over-damped dynamics and, second 
and most importantly, to prove the realizability of ([2]) and 
©• This means to determine the protocols for which ([2|) , 
(|3|) are attained. A theoretical effort in this direction was 
explicitly demanded by the authors of the micro-scale 
heat engine 

In the present contribution we address both questions. 
We specify in formula (1321) below the value of the con¬ 
stant a appearing in ([3]) in terms of the mobility for any 
reasonable probability distributions describing the state 
of the system at the end of the isothermal processes. 
Furthermore, we prove that power optimizing protocol is 
specified by the well-known Monge-Ampere-Kantorovich 
mass reconstruction algorithm @ . To neaten the no¬ 
tation, we report in details the calculations in the case of 
constant isotropic mobility, the extension to the general 
case being straightforward. 

We prove also a second universality result for the effi¬ 
ciency. It was noticed in 18| that the efficiency at maxi¬ 
mum power for a Stirling cycle well approximates also the 
efficiency at maximum work output in the limit of short 
duration of the cycle even for ideal heat cycles other than 
Stirling’s. Indeed, we show that if we optimize the work 
output with respect to the “expansion” state m,, and to 
the duration of the isothermal expansion then, in the 
limit of short cycle period (T i 0), the efficiency tends to 
o . We expect this result to be relevant for the design of 
nano-scale engines when controlling “target expansion” 
state nio may not be feasible. 

Let us briefly explain the origin of the universality of 
our results. It stems from the fact that the entropy 
production during an arbitrary thermodynamic transi¬ 
tion is proportional to the kinetic energy associated to 
the current velocity of the system. This fact was 
known for some time (see e.g. l^)- Its central role, 
however, in mapping optimal control of stochastic ther¬ 
modynamic transitions into optimal deterministic mass 
transport sm was only recently understood Q (see 


The system is driven by the gradient of a function 
?7 : R'’* X K I—>■ R which we take in the class of conhn- 
ing potentials sufficiently regular to justify the manipu¬ 
lations which follow. Let us also start by taking for a 
regularized, T-periodic and differentiable, version of the 
temperature time profile used in [^ . Accordin g to the 
general framework of stochastic thermodynamics |25Ll27l| . 
the mean heat dissipated by the working system during 
a temperature cycle is equal to the expectation value of 
the Stratonovich stochastic integral 

Q = -dh j\cfdc^U{Cut) ( 5 ) 

Well known properties of the Stratonovich integral yield 
after elementary manipulations 0 

E^%t|||u(C*,t)f+ 5(C,t)^^| (6) 

where v{z,t) = —dz{U{z,t)+l3^^S{z,t)) is the cur¬ 
rent velocity of the system, S{z,t) = — \n\^{z,t)/C] the 
microscopic entropy and St = dL S{C^,t). C is a dimen¬ 
sional constant irrelevant for the considerations which 
follow. From dH) we immediately see that if fit, = Pti, 
m(z, A) = m{z,tf) the only contribution to the dissipated 
heat comes from the entropy production. Owing to the 
non-degeneracy of the noise in the over-damped approx¬ 
imation, we are free to control the entropy production in 
terms of the current velocity v. Optimal control strate¬ 
gies are then most conveniently found by searching for 
extremals with respect to m and v of the Pontryagin- 
Bismut functional |17| 

A = E[V{Ctf,tt)-V{Ct,,Q] + 

Ej%i^\\viCt,t)r + - (SG)(Ct,t)|(7) 

The co-state function E: R^^ x R 1 —!> R plays here the role 
of a Lagrange multiplier imposing that 

{'DV){z,t) = [dt + v{z,t) ■ dz]V{z,t) (8) 

acts on scalars as a total derivative along the flow gener¬ 
ated by V. 
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Efficiency at maximum power. We construct the 
optimal cycle by optimizing m separately in , to) 
and [to,tf) under the boundary conditions = 

m{z,tf) = m,{z) = C exp{- 5 't( 2 :)} and 111 ( 2 :,to) = 
mo( 2 ;) = C exp {—S'o( 2 ;)}. We choose So and such 
that the probability densities are smooth, and that the 
Gibbs-Shannon entropy 

- St, = So-s, = E[^o(CtJ - > 0 (9) 

to signify an expansion of the system. In the limit of 
abrupt temperature changes, the derivative of the tem¬ 
perature vanishes except for Dirac-<5 contributions local¬ 
ized at to and tf. Proceeding as in , we find that 
along the isothermal branches of the cycle the optimal 
current velocity satisfies the two equations 

2)r; = 0 & TiS-d^-v = 0 (10) 

complemented by the stationarity condition dzV = 2v. 
The set of these three equations define the Monge- 
Ampere-Kantorovich optimal mass transport problem 
[^,[ 23 . Let us focus on the sub-interval [L,<o)- For any 
t € [ft, to) the velocity satisfies 


by m^{z) and 1110 ( 2 ;) [13|. Following [26| we define the heat 
input during the cycle as 

5o - 5. /C 

iT 

and the heat output 

So-S, 1C 


Qin — 


(16) 


Qout - 


Pc 


(i-7)r 
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The definitions imply —Q= Qm + Qout = Wo«t where 
Wout is the work output. The power of the cycle is then 


P = 


Wo 


So — s. 


K 


T 


= Vc- 


.r 


and its efficiency 


7(1 - 7 ) 


(18) 


_ Wout + (1-7) r 
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having used 1 — 7 c = Ph Pc^ ■ The maximum power is 
then attained for 7 * = 1/2 and 

SPhlC 
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v{cf){t;z,ti,),ti,) = v{z,t,) = 


^dzV{z,t,) 


( 11 ) 


with free steaming characteristics 


and it is equal to 

pI{So-S,)^ 
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( 21 ) 


cf){t-,z,t,) = z+ ^dzV{z,t,){t - t,) (12) 

The initial velocity is on its turn determined by the so¬ 
lution of the boundary problem 

S{z,to) 

= S{cj){to]z,t,),to) - tr\n[dz » (f){to; z,t,),to)] (13) 

In writing ()13p we imposed the continuity in time of the 
probability measure and of the Lagrangean map (1121) at 
to- The crucial observation is that since S{z,to) = S^z), 
S{(f>{to', z,ti,),to) = So{(p{to]z,tP}) are assigned, the La- 
grangian map 4){to\z,tP) = 4>i,{z) depends only on the 
boundary conditions and not on the duration of the 
isothermal process. As we can repeat the same consid¬ 
erations for the optimization in [to^f), we arrive after 
elementary manipulations (see e.g. [l| for details) at the 
general expression of the heat dissipated over one cycle 

For convenience we introduced the ratio 7 = {to—tP)/T € 
[ 0 , 1 ] and defined 

/C = E||c/,*(CJ-Cqf>0 (15) 


Finally evaluating the efficiency at maximum power 7 = 
7 *, T = 7/ yields ([ 2 ]) independently of nit and mg. The 
protocol attaining the maximum power is specified by 
which can be numerically solved using the algo¬ 
rithms given in [^, @ . This is the first of our announced 
universality results. 

Efficiency at maximum work output. We now con¬ 
sider a different optimization setting. As before, we 
allow the temperature to change only at to and tf. 
Also as before we assign m(z,tj = 111 ( 2 ;, tf) = m^z) = 
C exp {—£' 1 ( 2 :)} and we wish to optimize with respect to 
the duration of the isothermal processes. The difference 
is now that we look for minima of the dissipated heat 
(equivalently, maxima of the work output) with respect 
to m( 2 ;, to). 

In such a case, the optimal current velocity still obeys 
OT except when the isochoric temperature changes oc¬ 
cur. Upon requiring the continuity of the probability 
density at to we find 

v{z, to) - v{z, to) = vc 

We can then avail us of the continuity of the current 
velocity characteristics at to, their free streaming form 
along isothermal processes, and of (HU to derive the ex¬ 
tremal condition for the value of the initial velocity 


From the mathematical slant, K, is the squared Wasser- 
stein distance between the probability measures specified 


v{z,tP) 


(1 - 7 ) PC 


{dy,S){(f>{to;z,t,),to) 

2ph 


(23) 
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Combining this equation with m we finally get into 
a closed equation for the Lagrangian map (p{z) = 
z, ft) governing the maximum work output 

—T--r- ^ - dz tr ln(a 2 ® if) = dzS, (24) 

7(1-7)»7C / 


We complement this equation with the boundary condi¬ 
tion that Lp maps the support of the initial density into 
itself We emphasize that the stationarity 

condition 2v = dzV implies that the Lagrangian map ip 
is itself a gradient map. Hence, under suitable regularity 
hypotheses, we have reason to expect that the problem 
(1^ is well posed [^. We also observe that d^cp = 0 
generically solves the variational equation obtained by 
differentiating with respect to 7 if 7 = 1/2. In other 
words, solutions of (IMl) attain stationarity with respect 
to to if the isothermal expansion takes half of the cycle 
period. This result is intuitive in light of the fact that 
we could optimize the work output by first looking for 
the optimal duration of the expansion for fixed target 
state at to and then search for the optimal state. Finally, 
differentiating (l24l) with respect to T indicates that we 
cannot generically expect optimization with respect to T 
to be possible. We therefore conclude that the solution 
of (12411 describing maximum work output is obtained 
at 7 = 1/2 and depends parametrically upon T. 

We now are interested in solving (1241) in the limit of 
very short cycle period. As (HU holds up to to, we obtain 
for T and z sufficiently small 

P^-z + ^^^0-dzS, (25) 


If the initial value of the microscopic entropy confines 
fluctuations in the bulk region where (1^51) holds, we are 
then in the position to compute the efficiency of the cycle 
in short period limit. Namely, straightforward manipu¬ 
lations show that also in the present case the heat input 
and output are respectively amenable to the form (fTBl) 
and dni). Using (1251) we obtain within accuracy the ex¬ 
pression of the Shannon-Gibbs entropy 

Etrln (CtJ 

(26) 

° Ph JWLd oPh 


Similarly, we obtain for the squared Wasserstein distance 
the estimate 


K. 





(27) 


Upon inserting these results into the expression of the 
efficiency we obtain 


limry* 

rto 


2r]c 
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(28) 


We have therefore validated in the framework of over¬ 
damped stochastic thermodynamics the result of [l^ . 


Accordingly, the efficiency of work at maximum output 
can be approximated for small cycle period and small 
temperature ratio by same formula describing efficiency 
at maximum power. 

Elementary Extensions. In principle, anisotropic ef¬ 
fects in the solvent may require to consider 


dCt =-Mt • 5^^C/(Ct,t)dt-b y • dwt (29) 

Here Mt is a strictly positive definite symmetric ten¬ 
sor modeling the mobility of the solvent. For many ex¬ 
perimental applications, it is adequate to assume that 
the mobility tensor is constant in space and depends on 
time only because of changes in the solvent temperature: 
Mt = M(/3t). Under these hypotheses, we can repeat 
step-by-step the above calculations. The relevant indica¬ 
tors along isothermal transformations become the square 
Wasserstein distances : 


1C, = E {(0*(CtJ - CtJ • • (0*(CtJ - CtJ} (30) 

where = M(/3^), £ = h,c and is either the solution 
of the Monge-Ampere-Kantorovich equation or of the 
counter-part of ([24]) : 


2 Ph {dz G p) _ ^ ^ p- z 

7707(1-7) T 


dz tr \n{dz ® p) 


dzS (31) 


for C = 7 Mg ^ -I- (1 — 7 ) M,j In both cases, we obtain a 
generalization of ([ 2 ]) and (|28)) amenable to the form © 
with 


0 < 




y/iCh -b 


< 1 


(32) 


Using the lower and upper bounds for a, we readiliy re¬ 
cover the inequalities 


VC . , 


Vc 

2.-VC 


(33) 


derived in 11 | by assuming a time asymptotic behav¬ 
ior of the entropy production along isothermal processes 
inversely proportional to their duration. This is not sur¬ 
prising because the entropy production by optimal pro¬ 
tocols in over-damped thermodynamics exactly satisfies 
such condition [l|, [l^ . 


Conclusions and Perspectives. To summarize, we 
derived the explicit expression ([3]) and (I33|) of the 
efficiency at maximum power for the Stirling heat 
cycle in the framework of over-damped (Langevin- 
Smoluchowski) dynamics. We also determined the al¬ 
gorithm for computing the protocols achieving such effi¬ 
ciency. Present day technological advances occur at an 
impressive speed. Laboratory evidences of the possibil¬ 
ity of implementing Stirling cycles at the nano-sc&\e is 
already documented in the literature [l 6 | and the de¬ 
sign of a possible experiment has been recently proposed 
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0 , see also [3, [l^- These advances pose the chal¬ 
lenge to repeat the present analysis in the under-damped 
(Langevin-Kramers) regime. We recently proposed a 
general theory of optimal control of the dissipated heat 
for nano-mechanic al sy stems governed by a Langevin- 
Kramers dynamics 


22|. There we showed that the over¬ 


damped approximation is not only useful for perturba¬ 
tive analysis of the Langevin-Kramers but also provides 
a priori lower bounds on the dissipated heat. Based on 
these results, we expect that the present contribution will 
be relevant also for the analysis of nano-scale cyclic heat 
engines. 
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